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Detecting Hidden Differences via Permutation Symmetries 
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The predominant paradigm of quantum information science 
is the qubit, a quantum two-level system. This useful notion 
has allowed many important concepts to be abstracted away 
from particular physical implementations, revealing an under- 
lying structure in the way that information is manipulated and 
measured in quantum mechanics. Qubits are usually realised 
in some degree of freedom of a physical system. In many 
systems such as trapped ions and nuclear spins, the physical 
particles are inherently separated and the quantum statistical 
nature of the particles, be they bosons or fermions, can safely 
be neglected. In other systems including quantum optics, de- 
generage atomic gases, optically trapped atoms, and quasi- 
particles such as polaritons, the quantum statistics of the par- 
ticles can often play a role in the system's behaviour. Some- 
times the quantum statistical behaviour can be very useful, as 
in the Hong-Ou-Mandelfll effect in quantum optical systems 
which is often used to post-selectively implement interactions 
between photons JH. 

If a quantum information experiment is set up so that each 
particle is uniquely different in some observable degree of 
freedom - photons occupying different arms of an interferom- 
eter, say, or ions in different locations in a trap - then the quan- 
tum statistical properties of the particles generally do not play 
a role. Characterization of the quantum state of such systems 
proceeds according to the well-known procedures of quantum 
state tomography Ja|. The influence of external, unobserved 
degrees of freedom can be accounted for in this character- 
ization and results in a density matrix displaying less-than- 
perfect coherence. 

Recently, there have been several proposals and experi- 
ments involving multiple photons occupying a single spatio- 
temporal mode and two polarization modes [0,0, 01 ■ Such sys- 
tems do not fit into the qubit paradigm and quantum statistics 
plays an integral part in their behaviour. The states of these 
systems are of enormous interest ll27ll because they have been 
shown to exhibit phase superresolution in interferometry[0 
8], to be capable of beating the diffraction limit in lithog- 
raphy ji, , IToh . and to open up new avenues in quantum 
imaging ltlU U2I1 . They have also been proposed as a conve- 



nient qutrit useful in certain quantum cryptography and quan- 
tum information applications 10, HH. While, to date, pho- 
ton systems are the only ones where such states can be cre- 
ated, recent developments in optical lattices lflitl and else- 
where promise to open up similar opportunities in other phys- 
ical systems in the near future. 

Because these states involve multiple occupancy of a single 
mode, the quantum statistical nature of the particles is cru- 
cial to understanding their behaviour. Usually in considering 
such states the formalism of creation and annihilation opera- 
tors on the field mode is used. For example, the NOON state 
\N,0 : 0,N)OI can be written as 



10} , 



(1) 



where the subscript indexes the distinct modes. When such 
states are created experimentally, a central task is to recon- 
struct a faithful characterization of the state from measure- 
ment statistics. 

Ideally one would prefer to assume nothing about the 
source of the quantum states, treating it as a 'black box', and 
assume only that one has a set of measurements that one is 
able to accurately perform on a particular degree of freedom 
such as polarization. The reconstruction of the state from the 
measurements is called quantum state tomography and it has 
been an essential tool in quantum state engineering, quantum 
information science and quantum computing If the source 
produces an indefinite number of photons then continuous 
variable homodyne tomography methods can be extended to 
these states (1 5] . If the number of photons is known, though, it 
is simpler to extend the quantum state tomography techniques 
developed for qubits to systems of multiply occupied modes, 
as was done for example by Bogdanov et a/ lfl6ll . 

In their procedure one creates a basis of states from creation 
operators for a single spatio-temporal mode and the polariza- 
tion modes that the state can occupy. For the two-photon case 

that they studied their basis states were ^a H a' H ,a* H a\, ,aj,ay j , 

all taken to act on the vacuum. 
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There is a subtle assumption in writing the states this way 
that runs counter to the purpose of quantum state tomography. 
The notation for the states assumes that the raising operators 
act on the same spatio-temporal mode, but this is not some- 
thing that can be experimentally verified from the polarization 
measurements performed in the tomography. Indeed since the 
experiment in iflal involved combining different spontaneous 
parametric downconversion sources into a single spatial mode 
there is every reason to think that the different phase-matching 
conditions in the two downconversion crystals would result in 
the raising operators for each source acting on a somewhat 
different frequency-time mode. 

If the two raising operators do indeed act on different 
spatio-temporal modes then there can be a direct impact on 
polarization measurements since the amplitudes for differ- 
ent polarizations will carry the bosonic enhancement factors 
that one obtains by multiplying raising operators on the same 
mode. In principle one could attempt to fully characterize the 
spatial and temporal degrees of freedom to obtain the correct 
raising operator for each photon. Such a full characteriza- 
tion is technically very difficult, if not impossible. Moreover, 
full information about the spatio-temporal modes is likely not 
even desireable when it is ultimately polarization that is the 
degree of freedom of interest. 

What one would like is a 'black box' tomography technique 
for reconstructing the state in terms of polarization measure- 
ments only. The resulting description of the state ought to 
somehow include the influence of all unobserved degrees of 
freedom on polarization measurements. It ought also to cor- 
rectly predict the outcome of any polarization measurement 
one wishes to perform so as to be considered 'tomographically 
complete' . This paper develops and analyzes exactly such a 
technique. 

To our knowledge the problem of charactizing the state of 
multiply occupied modes has only arisen experimentally in 
photonic polarization systems 12811 . While we will concentrate 
on this specific realization, our method is completely general 
and can be applied to any of the aforementioned physical sys- 
tems in which quantum statistics play a role, either bosonic 
or fermionic. In order to have the discussion that follows re- 
flect this generality we will define some technical language. 
We will call the information-carrying degree of freedom in 
such systems the 'visible' degree of freedom. In the case of 
photonic polarization systems such as in JH [HI] this is polar- 
ization. All other degrees of freedom to which the apparatus 
is not sensitive we call 'hidden'. The description of the state 
that results from our state tomography procedure we label the 
'accessible' density matrix p acc . 

While inspired by practical problems encountered in our 
attempts to characterize quantum states, our approach is inter- 
esting in its own right as an exploration of how adding distin- 
guishing information in experimentally inaccessible degrees 
of freedom affects the quantum statistical properties of states. 

The key to our approach is to separate the state explicitly 
into hidden and visible parts and to examine the constraints 
placed on the visible degree of freedom by the quantum sta- 
tistical requirements on the whole state. For photons, the 
bosonic statistics require that the whole state be invariant un- 



der exchange of all particle labels. The exchange symmetry 
is more easily studied in state notation rather than raising op- 
erator notation, so we will use state notation throughout this 
paper. This can be confusing at first because often in the lit- 
erature states are written in a way that does not make the ex- 
change symmetry explicit. For example, one might write the 
polarization state of two photons as \HV), which is not obvi- 
ously exchange-symmetric as it must be for bosons. In such a 
description the order of the two labels implies the existence 
of a degree of freedom other than polarization, say differ- 
ent spatial modes a and b. a and b could be, for example, 
the distinguishable output angles of downconverted photons 
in spontaneous parametric downconversion. The full bosonic 
state needs to be symmetric under exchange of both spatial 
and polarization degrees of freedom, and would be written as 
| Vi) = (\HV) \ab) + \VH) \ba)) /\/2. Note that exchange of 
both the spatial and polarization labels leaves the state invari- 
ant, but the state has the property that one of the spatial modes, 
a is always correlated with one polarization H and the other 
mode b is always correlated with polarization V. In cases 
where the individual photon polarizations can be treated as 
qubits because a and b are distinguishable paths, this notation 
is redundant because no use is made of the permutation prop- 
erties of the whole state. For this reason it is usually preferable 
to write the state as \H a Vb) which denotes the correlation be- 
tween spatial and polarization modes without making explicit 
the bosonic exchange symmetry. It should be understood that 
in all circumstances this way of writing the state is simply a 
shorthand for \ ^\). 

We emphasize this point because we wish to discuss situa- 
tions where the overall exchange symmetry of the state is im- 
portant and the notation of \^f\) becomes very useful. There 
are situations where the spatial modes in the above example 
are 'hidden', that is to say they are not resolved by the detec- 
tion apparatus. This might occur if the photons were nearly 
collinear, but with a small angle between them. A multimode 
collection system such as a lens focusing onto a photodetec- 
tor significantly larger than the optical wavelength would have 
no means of distinguishing these two slightly different spatial 
modes. More generally, there could also be unresolvable 'hid- 
den' time-frequency modes that can become occupied due to 
uncorrected delays or dispersion. Since the nanosecond-scale 
resolution of most single-photon detectors is much longer than 
the femtosecond timescale of pulsed experiments, different 
time-frequency modes are generally not resolved by detectors. 
In non-photonic systems there are also myriad reasons why a 
given degree of freedom might be 'hidden' from experimental 
measurements. When a hidden degree of freedom is differ- 
ent for two particles we sometimes say that the particles are 
distinguishable in principle but not in practice. 

We can express our ignorance about the state of these 'hid- 
den' degrees of freedom by tracing over them. This leaves a 
density matrix observable only in the visible degrees of free- 
dom that we call the accessible density matrix: 

Pace = Tr hid [p] . (2) 

For example, if in the state 1 y\ ) the modes a and b cannot 
be resolved then we trace over them to obtain the accessible 
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density matrix 

pacc=Tr hld [|y/i)(V'i|] (3) 
= 1 -\HV)(HV\+ 1 -\VH)(VH\ (4) 

This is a mixed state of polarization. If the two photons 
had occupied the same spatial mode so that the state was 
(\HV) \aa) + \VH) \aa)) j \[2, then tracing over the spatial de- 
gree of freedom would have yielded a pure accesible density 
matrix in polarization \ (\HV) + \VH)) ((HV\ + (VH\). Since 
these two situations yield different density matrices on the 
polarization degree of freedom they can be distinguished by 
polarization measurements alone. The particular feature that 
distinguishes them is the antisymmetric part, expressed as the 
population of the singlet state (\HV) — \VH}) /VZ The sin- 
glet state projection makes up one element of the accessible 
density matrix. It is a measurable quantity even when the ex- 
perimental apparatus cannot tell the two photons apart. As 
we discuss extensively in lfl7tl . for two photons the presence of 
an antisymmetric component of the polarization state implies 
the existence of one or more unobserved degrees of freedom 
that are different for the two particles and correlated in just the 
right way to result in the correct bosonic exchange symmetry 
for the whole state. This shows that differences in the hid- 
den degrees of freedom may be inferred from measurements 
performed only on the visible degrees of freedom. 

The remainder of this paper will examine how the acces- 
sible density matrix can be calculated and measured for an 
arbitrary number of particles and for a visible degree of free- 
dom with an arbitrary, finite number of levels. In the first 
section we will begin by determining how many elements are 
contained in a general accessible density matrix as a function 
of the dimensionality of the visible degree of freedom and the 
number of particles. This determines both how many linearly 
independent measurements can be made and how many num- 
bers are needed to calculate all possible expectation values on 
the visible degrees of freedom. The second section will put the 
discussion of the first section on a firm group-theoretical foot- 
ing. In the third section we examine how the theory applies to 
the case of three photon polarizations. In section 4 we discuss 
how the accessible density matrix can be measured and work 
through a specific numerical example with three photons. Fi- 
nally, in section 5 we discuss what claims can be made about 
the indistinguishability of the particles from a knowledge of 
the accessible density matrix. 

I. THE FORM OF THE ACCESSIBLE DENSITY MATRIX 

In this section we develop the structure of the accessible 
density matrix and show how many independent measure- 
ments can be done on a visible degree of freedom. We will 
start by assuming a two-level degree of freedom like photon 
polarization and then extend the result to a of-level degree of 
freedom that could, for example, be the Laguerre-Gauss spa- 
tial modeJTH of photons. 

Our approach is to consider the Hilbert space of the photons 
as a tensor product of a Hilbert space describing the visible de- 



grees of freedom and another describing the hidden degrees of 
freedom. Consider Af photon polarizations. Polarization trans- 
formations e I { xa "+y a y +za -~) (where {o x ,O y ,O z } are the Pauli 
matrices) give a realization of the group SU(2). Since SU(2) 
acts irreducibly on the photon polarization, we can view the 
photons as spin one-half systems; that is, they transform ac- 
cording to the 7 = 5 irrep of SU(2). If all systems are dis- 
tinguishable in practice, so that each photon is in a separate 
mode that can be experimentally distinguished (different rails 
of a multi-rail interferometer, say), then the dimension of the 
accesible space is the full dimension 2 N and the number of 
accessible density matrix elements is 2 2N . This is the familiar 
situation of quantum state tomography as applied to photon 
polarization 11 811 , trapped ions and other qubit systems. We 
would like to know the comparable number of density ma- 
trix elements when the photons are not experimentally dis- 
tinguishable because the degrees of freedom that might dis- 
tinguish the particles cannot be resolved experimentally. We 
note that in this case 2 2N provides an upper bound on the num- 
ber of elements in the accesible density matrix. 

The following decomposition of the Af -polarization Hilbert 
space will be useful. Unitary polarization operations acting on 
the whole state can be decomposed into 'angular momenta' j 
(irreducible representations of SU(2)) according to the well 
known Clebsch-Gordan series; for example 

- 2 = ieo, 



- = 2ffilffil©lffi0©0. 

Notice that if N > 2, certain j values occur more than once; 
they are said to have multiplicity. However the largest j al- 
ways occurs only once, since there is only one way to couple 
the spin- A particles to maximum j = The states in the ^ 
space are always totally symmetric under permutation of the 
N polarizations. If they are indistinguishable in principle, i.e. 
their hidden degrees of freedom are in the same state, then 
these totally symmetric visible states are the only ones avail- 
able to the whole state by the restriction that it have bosonic 
symmetry. Since the dimension of a spin j space is 2j + 1, 
in this case the dimension is 2^ + 1 = N + I and the num- 
ber of accessible density matrix elements is (N +i) 2 . Pre- 
vious tomography schemes such as the one used in 10, fl6tl 
worked under the tacit assumption that the photons were in- 
distinguishable in principle, and so described the polarization 
only in terms of these j =N/2 states. 

For experimentally distinguishable particles we see that the 
number of density matrix elements grows exponentially as 
2 2N with the number of particles. And as we have just shown, 
for particles indistinguishable in principle, the number of ele- 
ments grows polynomially as (N +1) 2 in the number of par- 
ticles. How does the number grow when the particles are dis- 
tinguishable in principle, but not in practice? In this case, we 
must trace out the hidden degree of freedom in order to ex- 
press our ignorance about them, but in doing so we are forced 
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to erase the ordering information of the N systems. This infor- 
mation is encoded both in the phase between different terms 
in the Clebsch-Gordan decomposition and, when multiplicity 
is greater than one, in how population is distributed among the 
orthogonal eigenvectors of the multiplicity space; in terms of 
operations, the unitary polarization transformations take states 
with angular momentum j to other states with angular mo- 
mentum j in the same multiplicity space, while permutations 
take states from one multiplicity space to a different multiplic- 
ity space of the same j. 

Sectors of states all carrying the same value of j form 
(2j + 1) by (2j + 1) block-diagonal submatrices along the 
main diagonal of p acc . SU(2) operations rotate states within 
these blocks and permutations of the polarization labels move 
population from one block to another with the same value of 

j- 



* * * * 
$ $ * * 

* $ * * 

* * * * 
Sat acts 

between /* 
blocks \ 



SU(2)acts 

within 

blocks 



(6) 



This explains why the highest j space is symmetric — the 
space has multiplicity one, and so must be invariant under per- 
mutations. 

When we trace out the hidden degrees of freedom, coher- 
ences between states of different j as well as all information 
about the state within the multiplicity spaces are destroyed, 
leading to a density matrix that is block diagonal. A conse- 
quence is that populations in multiple copies of the same j are 
averaged, yielding multiple copies of the same density sub- 
matrix. Thus the accessible density matrix consists of only 
one independent density submatrix for each j in the Clebsch- 
Gordan decomposition with zero coherence between subma- 
trices. The number of independent accessible density matrix 
elements is therefore 



I (2i- 



j=0 or 



N + 3 



(7) 



Thus the number of density matrix elements scales polynomi- 
ally in the number of particles, at least for two-level systems 
like polarization. 

When the visible degree of freedom has d distinct levels the 
situation is completely analogous, with the Clebsch-Gordan 
series generalised to SU(d). The space of N <i-level sys- 
tems decomposes into irreps A of SU(c/) where now the label 
A = {X\,Xz,-- - ,Xa); = N, A; > A,-+i is a regular par- 
tition of N (a Young diagram). Of course, if the systems are 
experimentally distinguishable then the entire d N dimensional 
Hilbert space is accessible and the number of accessible den- 
sity matrix elements is d 2N which gives the (exponential) up- 
per bound. The irrep (N,Q,0, ■ ■ ■ ,0), (analogous to highest j 



in the SU(2) case), occurs only once in the decomposition of 
the Hilbert space and so is always symmetric under permuta- 
tions. The dimension of this (and indeed any SU(t/)) irrep is 
given by the Weyl character formula lfl9tl 



dim(A) = 

\<i<j<d 



Xj-Xj + j -l 
j-i 



(8) 



If the qudits are indistinguishable in principle, then again they 
are restricted to the totally symmetric subspace with Ai = N 
and all other Xj = 0; the Weyl formula gives 



N - 



n 

j=2 J 



J 



1 



N+d- 
N 



(9) 



for the dimension, so the number of accessible density matrix 

elements is f + ^ 1 ) 2 - 

The unitary and permutation group actions are the same as 
in the SU(2) case. SU(cf) acts within irrep spaces A and Sat 
acts 'across' multiplicities. When the distinguishing degrees 
of freedom are hidden, the ordering information of the sys- 
tems is lost and the permutation group action is trivialised, 
leaving only one 'copy' of each SU(c/) irrep space for each A . 
The dimension of the accessible space is therefore |29|| 



i n 

X \<i<j<d 



Xj - Xj + j - 



N + d 2 -I 
N 



(10) 



and is always a polynomial in A^. 



II. GROUP THEORETICAL CONSTRUCTION 

Here we will construct explicitly the most general totally 
symmetric state of a system of particles with both visible and 
hidden degrees of freedom, which we use to justify the claims 
made in the last section. Let the Hilbert spaces for these two 
be denoted Jt° vls and J$? hld , respectively. Assuming that there 
are Af particles, the same permutation group Sat acts on both of 
these spaces. Decompose each space into irreps A of Sn and 
consider their tensor product, the space of all available states: 



Jif. 



hid 



(ID 



, I'. i 



where m labels the multiplicity of irrep A. 
orthonormal basis for each irrep space J%\,r, 
are labelled 



|Am M ) vis |A'mV'; 



hid ' 



Let jl index an 
; the basis states 



(12) 



where m = 1,2, ■ • • , mult^^ runs over the multiplicity of irrep 
A in the Hilbert space, and jl = 1 ,2, • • • , dim,J£| runs over the 
dimension of irrep A . For readers familiar with Schur-Weyl 
duality, m indexes a basis for an irrep A of the unitary group 
action on each particle, and jj. indexes a basis for an irrep A 
of the permutation group action Sn- The fact that the same 
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irrep label can be used for both group actions is why they are 
'dual'. 

Now the problem of finding totally symmetric states in ^ 
is a coupling problem, completely analogous to coupling an- 
gular momentum states to arrive at states of angular momen- 
tum zero. In fact, it can be shown from the rules for tensor 
products of Young diagrams that the totally symmetric irrep 
A = (N) of Sjv only occurs in a tensor product A (g) A' if A' = A, 
and moreover that (N) only occurs once, i.e. it has multiplicity 
one l2(ill . The analogy is that the spin zero irrep of the rotation 
group only occurs in the tensor product j £g> / if / = j, and it 
occurs only once, i.e in order to couple two angular momenta 
to j = 0, we know the two angular momenta must be equal. 
Note also that (N) is always one dimensional, so there is one 
totally symmetric state for each A, unique up to multiplicity. 

Given an irrep A and two multiplicity sectors m,m' in Jrff, 
this unique totally symmetric (unnormalised) state \Xmm') is 
an equally weighted superpositon of the states of each factor 
in the tensor product 

\Xmm') = £ |Am^) vis |Am' i u}hid> ( 13 ) 
n=i 

(which is a state on the combined space, not to be confused 
with the uncombined visible and hidden states, despite the 
fact that they both have three labels). The most general to- 
tally symmetric pure state in Jtf 1 is therefore an arbitrary linear 
combination of these: 

\Yn)=ZL C L'\^')- (14) 

X mm! 

The same analysis goes through for totally antisymmetric 
states. The unique coupling is A ®X, where A is the irrep 
conjugate to A. There is a restriction, however, given by 
the dimension of the Hilbert space for each particle, which is 
again encoded in the rules for Young diagrams. For example, 
there is no totally antisymmetric state of three indistinguish- 
able spins. 

Now we can define what we mean by distinguishable and 
indistinguishable. Expand | y/jy) in the physical basis of N par- 
ticles. Those states in the expansion where the hidden degrees 
of freedom for all N particles are in the same state are indistin- 
guishable in principle. This hidden state is totally symmetric 
by definition, and by the coupling mentioned above it follows 
that the visible state must also be symmetric. Since (N) is one 
dimensional, there is only one term in the sum over the basis 
indexed by jj. above, and the total state is separable across the 
hidden and visible subspaces. Thus, tracing out the hidden 
space does not alter the visible state, and since it can only lie 
in (N), the accessible density matrix is restricted to the totally 
symmetric subspace, as expected. 

Those states in the expansion where the hidden degrees of 
freedom for all N particles are in distinct orthogonal states 
are distinguishable in principle. There is a large amount of 
entanglement across the hidden and visible subspaces. If the 
hidden modes are inaccessible in practice, then we arrive at 
the accessible density matrix by tracing out the hidden modes. 



Using Eq.(fT3ll. one finds 

Pa C c = Tr hid [|v^)(v^v|] (15) 
= £(cr£v| £ C^C^lXmm^iKnn'Walv^ 

Gtv X Kmm'nn' 

(16) 

= I PwII^M^/nVl. (17) 

Xmm' M 

One therefore concludes that p x , = V„Ci„C^f affords the 

"mm ^ n mn mn 

only freedom in the accessible density matrix, giving only 
one value per irrep A and pair of multiplicity indices m,m'. 
The trace erases coherences between different A sectors on 
account of those sectors being orthogonal. We also see that 
the equally weighted average over jj. which was necessary for 
total symmetry has destroyed any independence between the 
multiplicity spaces — we get the same copy of the A subma- 
trix for all fi, and so we effectively have one submatrix for 
each A. From the point of view of accessible measurements 
the state space has 'collapsed', although if the particles were 
distinguishable in one of the hidden degrees of freedom, then 
the ability to measure that degree of freedom would restore 
the Hilbert space to its full size. 

The measurement of the accessible density matrix elements 
C^ m allows one to infer the existence of hidden differences 
among the particles making up the state. To see this, consider 
that the hidden and visible spaces must both transform under 
the same permutation group S#. If we decompose the visible 
and hidden spaces separately under this common group action, 
we arrive at visible states labelled by Sn irreps A and hidden 
states labelled by Sat irreps A'. Again, coupling visible and 
hidden states to make totally symmetric states is completely 
analogous to coupling angular momentum states to make an- 
gular momentum j = 0. It follows that A' must equal A. 
Thus, if a visible state is measured to be in a state of per- 
mutation symmetry A that is not totally symmetric, one can 
infer that there existed a hidden state of permutation symme- 
try A to which it was coupled, implying in turn the presence of 
multiple orthogonal states for the hidden degrees of freedom. 
These hidden differences serve to make the the photons distin- 
guishable and explain why the coherences between different 
A (j for SU(2)) disappear when the hidden states are traced, 
simply because states of different A are orthogonal. 



III. EXAMPLE: THE ACCESSIBLE DENSITY MATRIX 
FOR THREE PHOTON POLARIZATIONS 

To make the discussion of the previous sections more con- 
crete we will focus on the particular example of three pho- 
ton polarizations. This example is experimentally relevant to 
previously published work from our group on NOON states [U, 
and to ongoing work on making other states in the same three- 
photon polarization Hilbert space. 

The Clebsch-Gordan decomposition for three spin-^ parti- 
cles was given in Eq. ([5j. We can explicitly write out the 
states of this decomposition. Each state is labeled by a pair of 
angular momentum quantum numbers j and m. The j = 3/2 
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FIG. 1 : Apparatus for measuring the accessible density matrix for N 
photon polarizations. The state I 1 ?) is sent into a quarter waveplate 
(QWP) and half waveplate (HWP) followed by a polarizing beam- 
splitter (PBS). Number resolving photon counters count the number 
of vertical photons Ny and the number of horizontal photons Nh- 

states that are completely symmetric under permutations are: 

|3/2, 3/2) = \HHH) , (18) 

\/3|3/2,l/2) = \HHV) + \HVH) + \VHH) , (19) 

V3\3/2, -1/2) = \VVH) + \VHV) + \HVV) , (20) 
|3/2, -3/2) = \VVV) . (21) 

While the j = 1/2 space has multiplicity two. The two spaces 
are spanned by 

V6 1 1/2, 1/2) 1 = \HHV) + \HVH) - 2 \VHH) , (22) 

V6 1 1/2, -1/2)! = \VVH) + \VHV}-2\HVV) , (23) 

and 

V2 1 1 /2 , 1 /2) 2 = \HHV ) - \HVH) , (24) 

V2\l/2,-l/2) 2 = \VVH}-\VHV}. (25) 

1 1/2, 1/2) j transforms into 1 1 /2 , 1 /2) 2 under permutation op- 
erations and 1 1 / 2 , — 1 J 2) l transforms into 1 1 /2 , — 1 /2) 2 in ex- 
actly the same way. However, polarization measurements can- 
not distinguish 11/2,1/2)! from 1 1/2, 1/2) 2 or 1 1/2, -1/2) l 
from 1 1/2, — 1/2) 2 . All they can do is determine the average 
of the two-by-two density matrix over the space spanned by 
1 1/2, l/2) 2 anc ^ 1 1/2, l/2)i anc l tne density matrix over the 
space spanned by |l/2, — 1/2) 2 and |l/2, —1/2),. From the 
point of view of polarization measurements, the information 
contained in the two spaces collapses into a single effective 
j = 1/2 sector of p acc . 

The accessible density matrix contains ( 3 ^ 3 ) = 4 2 + 2 2 = 
20 elements. When distinguishing information is hidden, the 
best characterization of the state of three photon polarizations 
is the determination of these 20 elements. 

IV. MEASURING THE ACCESSIBLE DENSITY MATRIX 

We have shown that elements of the accessible density ma- 
trix offer the most complete description of the state of N par- 
ticles when one degree of freedom of the particles is visible 



and others are hidden. It is not clear from our discussion so 
far that it is possible to measure p acc using available experi- 
mental tools. In this section we will show that in the case of 
polarization it is indeed possible to measure p acc with a sim- 
ple experimental device. This device, shown in figure Q] in- 
volves four different optical elements, a quarter wave-plate, a 
half-waveplate, a polarizing beamsplitter (PBS) and number- 
resolving photon counters (such as the one demonstrated in 
l2~lll ) at the two ports. The multiphoton polarization state 
passes through the two waveplates and on to the PBS where it 
is split into H and V components. The number of photons in 
each port is then measured with the number-resolving single 
photon detector. Such a device is completely analogous to the 
Stern-Gerlach device for measuring spin projections. 

If there are N photons in the state then there are N + 
1 different ways that these can split between the H and 
V ports of the polarizing beamsplitter. When the pho- 
tons leaving each port are counted, the measurement im- 
plemented will be a convex sum of projectors onto all 
states having that number of horizontal and vertical pho- 
tons. For example, P N = \HH - H) {HH ■ - H], P N -i = 
(\VH ■ ■ -H) (VH ■ ■ -H\ + \HV ■ ■ -H) (HV ■ ■ -H\ + . . .), and so 
on. In angular momentum language is a pure projector 
onto the state j = N/2, m = N/2. By changing the angles 
of the waveplates one can 'orbit' this measurement in the 
j = N/2 space thereby obtaining all the density matrix ele- 
ments lying in the j = N/2 space. All j — N/2 matrix ele- 
ments can be determined by measuring only rotated versions 
of P^. P;v_i is a convex sum of projectors onto all states with 
m = N/2—\. Since the j —N/2, m =N/2 — 1 projection can 
be determined from rotated versions of P#, we can subtract 
this part from Pn-i leaving a projector with support only in 
the j = N/2 — 1 subspace. By changing the waveplate angles 
one can use this reduced operator to completely characterize 
the j =N/2— 1 space. One can then subtract the j = N/2 and 
j = N /2— \ terms from the m=N/2 — 2 operator, and so on. 
In this way all the terms in the accessible density matrix can 
be measured. 

It should be noted that this is not the only way to measure 
the accessible density matrix. In fact any set of linearly in- 
dependent measurements equal in number to the number of 
accessible density matrix elements as calculated from equa- 
tion|7]can be used to reconstruct the accessible density matrix. 
Standard inversion techniques such as linear inversion, max- 
imum likelihood fitting flal and convex maximum likelihood 
fitting ly can be used to obtain p acc from an experimental 
dataset. 



A. Three-photon example 

In this section we work through an example of how our 
techniques might be applied in a three-photon polarization ex- 
periment. We consider an experimental situation similar to 
the one used to create the state -4- (|3#,0y) + |0#,3y)) in ref- 
erence 01. There three photons were combined on a beam- 
splitter and post-selection was used to isolate those instances 
where all three photons left from the same port of the beam- 
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splitter. If all the photons were indistinguishable, and each 
photon was set to the correct polarization then by this proce- 
dure the entangled state |3 , : 0, 3} = 4= (\3 H , 0y } + 1 H , 3y )) 
would have been. This comes about because of the state is 
factorizable in raising operators through the relation: 



( 4 + a V ) 



i/3 4) (4 



+ 



(26) 



Note that the right side is a product of polarization raising 
operators all acting on the same spatio-temporal mode. 

In that experiment, however, two of the photons were pro- 
duced by a spontaneous parametric downconversion process 
and the third was produced by an attenuated laser pulse. It 
is to be expected that these different sources might produce 
photons with hidden differences in their time-frequency wave- 
functions. In principle such differences can be reduced by fil- 
tering, but let us suppose that filtering is insufficient, resulting 
in the mode of the third photon having only a 50% overlap 
with the mode of the other two photons, which are identical 
to one another. We can model this by replacing the raising 
operators in the third bracket in equation [26] with operators 

4/y = 4^ ^4/y + 4/y) where is a creation operator for 

a mode b orthogonal to the mode of a for which is the 
raising operator. The 50% overlap is chosen here to keep the 
calculation simple. For a more general situation one can re- 
peat the analysis using 4/y = cos ^4/y + sin 0b r H ^ v where 
6 is an angle parametrizing the degree of overlap. 

Inserting this substitution into equation [26] we obtain the 
expression 

^(4+4) (4+^/34) (4+ ,4-734 



a H b H + a H a V b H I 1 



+ai 2 bie 2 «^ +4 2 4 e 4 *73 +444 f l+e 



Jix/3 



-OybX ] 



(27) 



Our goal is to write the state as an accessible density matrix, 
purely in terms of polarization measurements that can be done 
without knowing anything about the differences between the 
hidden time-frequency degrees of freedom. To arrive at such 
an expression we will need to trace over the orthogonal hidden 
modes a and b. To do so the expression must first be rewritten 
in a first-quantized notation that clarifies the imposed separa- 
tion between the hidden and visible degrees of freedom. The 
full expression is too long to write here, but rewriting one of 
the terms should be enough to give a feel for the calculation. 

Consider the term a* H a v b v . In rewriting this in first- 
quantized form we need to use tensor products of state vec- 
tors on the hidden and visible degrees of freedom, keeping in 
mind that the indistinguishability of the particles means that 
any one of the three can be in mode b. We write it as follows: 



444 



+ (\VHV) + 



+ \VHV)) \aab) + (\HVV) + \VVH)) | 
\VVH))\baa)} (28) 



The trace over the hidden degrees of freedom produces an 
incoherent sum over the density matrices for each bracketed 
polarization state since \baa), \aab) and \aba) are orthogonal. 
The resulting accessible density matrix describing this term is 



Pace — 








2/3 











1/6 





1/6 



(29) 



In the same way we can obtain the accessible density matrix 
for the entire state in Equationl27l 



■ 0.3636 


0.0 
0.3636 



0.3636 



0.3636 



0.0682 
-0.0341+0.0590/ 



-0.0341 -0.0590/ 
0.0682 



(30) 



0.0682 
-0.0341+0.0590/ 



-0.0341 -0.0590/ 
0.0682 



Note that the partial distinguishability of the third photon re- instead of the j = 3/2 spaces. The fidelity l23ll of this state to 
suits in 27% of the population being in the j = 1/2 spaces the desired state |3,0: 0,3) is 0.7273. The distinguishability 
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QWP HWP 3/2 1/2 -1/2 -3/2 



0° 


0° 


3645 1459 


1385 


3586 


15° 


0° 


2201 3953 


1006 


2703 


30° 


Q c 


275 7699 


160 


1932 


45° 


0° 


905 5260 


2904 


904 


C 


12.25° 


2078 2042 


3834 


1975 


15' J 


12.25° 


2759 2388 


2185 


2673 


30° 


12.25° 


2105 2693 


4174 


1108 


45- 


12.25° 


420 6700 


1459 


1272 


0° 


22.5° 


910 2741 


5163 


888 


15° 


22.5° 


892 4226 


3021 


1899 


30° 


22.5° 


1337 3838 


3207 


1550 


45° 


22.5° 


1914 2043 


6069 






TABLE I: Simulated results of measurement of the state 1 30 : 03) = 
4- (|3h,0v) + |0#,3y}) when the spatio-temporal wavepacket of 
one of the photons has a 50% overlap with that of the other two. 
The first two columns give angles for the waveplates, and the last 
four give the number of counts observed for each of the four out- 
comes of a number-resolving measurement. These outcomes are 
labeled by the value of m with the understanding that they include 
contributions from all spaces with j > m. For example the m = 1/2 
column corresponds to the measurement operator \HHV) {HHV\ + 
\HVH) (HVH\ + \VHH) (VHH\. 

of one of the photons can therefore make a significant differ- 
ence in the overall quality of the state. 

So far we have assumed that we know the exact behaviour 
of the hidden degrees of freedom for our state. Let us now 
instead assume the experimental situation in which we can do 
polarization measurements but do not know about the hidden 



degrees of freedom making one of the three photons different 
from the other two. We will use a detection apparatus like the 
one in figure [T] The quarter and half waveplates are set to the 
angles listed in the first two columns of table Q] This results 
in a number of detections for each measurement outcome as 
listed in the last four columns. The numbers were generated 
via Monte Carlo simulation of Poisson-distributed data arising 
from the density matrix in equation [30] It was assumed that 
on average 10,000 three-photon states were measured for each 
waveplate setting. 

The set of measurement operators in table U is overcom- 
plete, as can be verified by explicit calculations of the dimen- 
sion of the vector space they span. The 48 projectors span a 
space of 20 linearly independent dimensions. As predicted by 
equation[7] this is the maximum number of independent mea- 
surable operators when polarization is the only visible degree 
of freedom. 

These twenty parameters can be arranged to form an acces- 
sible density matrix in the form of equation|6] with the 20 ele- 
ments broken into a 16-element symmetric j = 3/2 subspace 
and the remaining four elements representing an average over 
the two j = 1/2 subspaces. 

Once this form is assumed for the accessible density matrix, 
the data can be fit to it using maximum-likelihood fitting lr22tl . 
To perform the fit we use the free convex optimization pack- 
age SeDumi lf24ll for Matlab. In order to measure the likeli- 
hood that a given density matrix gave rise to the dataset we 
calculate the log-likelihood l22ll . The density matrix that max- 
imizes this function given the outcomes listed in table U is 
given below. 



0.3626 0.0057+0.0033; 0.0001-0.0003; 0.3597+0.0010! 

0.0057-0.0033; 0.0036 -0.0006-0.0028; 0.0023-0.0040; 

0.0001+0.0003; -0.0006+0.0028; 0.0023 0.0013-0.0023; 

0.3597-0.0010; 0.0023+0.0040; 0.0013+0.0023; 0.3601 

0.0686 -0.0322-0.0597; 
-0.0322+0.0597; 0.0670 

0.0686 -0.0322-0.0597; 
-0.0322+0.0597; 0.0670 

(3D 



This can be seen to be very close the the density matrix in 
equation [30] with the difference accounted for by the statisti- 
cal noise in the measurements. 

The measured non-zero population in the non-symmetric 
subspace indicates the presence of hidden distinguishing in- 
formation. The detection of this population would allow an 
experimentalist to infer the presence of a hidden degree of 
freedom (in this case the time-frequency degree of freedom) 
distinguishing one of the photons from the other two. It will 
also be noted that to the extent that the photons are indistin- 
guishable, they are indeed in the desired state. In other words, 
the errors have arisen solely from the distinguishing informa- 
tion, and not, say, from unknown polarization rotations. This 



is all valuable information useful in diagnosing problems with 
the experiment. 



It should be emphasized that there is nothing special about 
the particular waveplate settings used in this example. The 
important thing is that the resulting measurement operators 
fully span the space of accessible density matrix elements. If 
this is the case then the maximum-likelihood problem is well- 
defined and guaranteed to converge to the unique solution Ir22ll . 
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V. USING THE ACCESSIBLE DENSITY MATRIX TO 
INFER FUNDAMENTAL DISTINGUISHABILITY 

One of the main reasons for characterizing an experimen- 
tally generated polarization state is to substantiate claims that 
a particular quantum state of light has been achieved. For 
nearly all quantum protocols only the j = N/2 symmetric 
states will be useful since all the other states involve unwanted 
correlations with the hidden degrees of freedom which, by 
definition, cannot be manipulated. Usually one makes the 
claim that all the photons in the state are 'indistinguishable' 
in the hidden degrees of freedom, meaning that they all oc- 
cupy the same hidden state. Our technique provides the first 
general method for verifying this claim. 

If, when the accessible density matrix is measured, all the 
population is found to be in the symmetric space then it must 
be true that the hidden degrees of freedom are also in symmet- 
ric states. If this were not true then the overall state could not 
have the requisite bosonic symmetry under permutation. 

If, in addition, the purity of the visible state is unity then the 
hidden degrees of freedom are unentangled with the visible 
degrees of freedom. This means that all measurements on the 
visible degrees of freedom will be consistent with the photons 
all being in the same single-particle hidden state. This being 
the case, it makes sense to call the photons 'indistinguishable' 
in the conventional sense of the word. 

This definition of indistinguishability is entirely consistent 
with the one proposed by Liu et al. |QJ, |25t] but is more flexible 
because it is expressed in the density matrix formalism which 
transforms in a predictable way under various operations one 
might wish to perform on the state. 

VI. CONCLUSIONS 

We have outlined a procedure for measuring the state of a 
system of particles spread over several experimental modes 



which may be entangled with hidden degrees of freedom. Our 
technique should be used to justify claims of production of 
'indistinguishable' photons. It is the most complete descrip- 
tion of the state possible when some degrees of freedom are 
hidden, and in particular it gives a more complete description 
of the state than previous characterization techniques such 
as those employed in fjfill . U] or ll25ll . In addition to being 
complete, this characterization also has the advantage of pro- 
ducing a density matrix that can be used in the usual way 
to predict the outcome of all measurements. We expect this 
method to become the standard means of characterizing states 
of a fixed number of experimentally indistinguishable photons 
just as quantum state tomography lfl8l has become the stan- 
dard means of characterizing distinguishable photons. Indeed 
since the number of accessible measurements for experimen- 
tally indistinguishable photons only grows polynomially with 
the number of photons in the state, our technique should prove 
useful for much larger systems of photons than state tomogra- 
phy does for distinguishable photons. 
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